Draft version February 27, 2013 

Preprint typeset using MTgX style emulateapj v. 5/2/1 1 



A REDLINE STARBURST: CO(2-l) OBSERVATIONS OF AN EDDINGTON-LIMITED GALAXY REVEAL STAR 

FORMATION AT ITS MOST EXTREME 

J. E. GEACH^, R. C. HICKOX^, a. M. DIAMOND-STANIC^, M. KRIPS*, J. MOUSTAKAS^, 

C. A. Tremonti", a. L. Coil^, p. H. Sell'', G. H. Rudnick^ 

Draft version February 27, 2013 



m 
T— I 

o 
in 

(N 



ABSTRACT 

We report observations of the CO(2-l) emission of SDSS J1506+54, a compact (re « 135 pc) starburst 
galaxy at z — 0.6. SDSS J1506+54 appears to be forming stars close to the limit allowed by stellar radiation 
pressure feedback models: the measured Lir/Lqq « 1500 is one of the highest measured for any galaxy. With 
its compact optical morphology but extended low surface brightness envelope, post-starburst spectral features, 
high infrared luminosity (Lir, > lO^^'^L©), low gas fraction (Mn^/Mi, ~ 15%), and short gas depletion time 
(tens of Myr), we speculate that this is a feedback-limited central starburst episode at the conclusion of a major 
merger. Taken as such, SDSS J1504+54 epitomizes the brief closing stage of a classic model of galaxy growth: 
we are witnessing a key component of spheroid formation during what we term a 'redline' starburst. 

Subject headings: galaxies: starburst, evolution 
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I. nvlTRODUCTION 

Theoretical studies suggest that the majority of the stellar 
mass in the central regions of massive spheroids forms in situ 
in powerful, compact starbursts (e.g. Narayanan et al. 2010). 
In these systems, radiative feedback from stars is expected to 
produce an upper limit on the star formation rate density (S^) 
by limiting the density of the star-forming gas (Sgas, Thomp- 
son, Quataert & Murray 2005; Murray, Quataert & Thompson 
2005). One explanation for the small range in central stellar 
surface densities observed in the cores of local spheroids is 
that they formed in compact starbursts where the bulk of the 
gas reservoir formed stars at the Eddington limit (Hopkins et 
al. 2010), however observations have only just begun to probe 

the most extreme S^-Sgas in distant powerful starbursts that 
are the antecedents of massive galaxies today. 

Diamond-Stanic et al. (2012, DS12) reported the discov- 
ery, with the Wide-field Infrared Survey Explorer (WISE), of 
a population of moderate-redshift (z w 0.5) galaxies that 
appear to be among the highest density starbursts yet ob- 
served. Hubble Space Telescope (HST) imaging reveals that 
these systems are remarkably compact in the optical bands, 
with re--100pc (DS12, R H. Sell et al. 2013 in prepa- 
ration), while the optical spectra and UV-IR spectral en- 
ergy distributions imply high star formation rates (SFRs) ap- 
proaching lO'^A/oyr^^ (and thus very high S*), with little 
or no contribution from an active galactic nucleus (AGN). 
The galaxies also exhibit extremely blueshifted (1000km s^^) 
Mg II AA2796, 2803 interstellar absorption lines (Tremonti et 
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al. 2007) indicative of powerful feedback. 

In order to test if these compact galaxies represent the for- 
mation of a massive stellar bulge via a feedback limited star- 
burst, we require a measurement of the cold gas reservoir, 
since Eddington limited star formation places a theoretical up- 
per limit on S^^/Sgas for the actively star-forming gas. In this 
Letter we present new observations at 2 mm of the most ex- 
treme galaxy in the DS12 sample with the Institut de Radioas- 
tronomie Millimetrique (IRAM) Plateau de Bure Interferome- 
ter (PdBI). We measure the CO(2-l) molecular emission line 
and derive the gas properties of this remarkable system to test 
the picture of Eddington limited star formation in a compact 
starburst galaxy. 

2. SDSSJ150636.30+540220.9 

SDSSJ150636.30H-540220.9 at z = 0.608 (hereafter 
'SDSS J1506-I-54') was targeted because it is the most ex- 
treme system in the DS12 sample, in terms of its star forma- 
tion rate density. The effective radius measured in HST 'WFC3 
F814W (rest-frame V-band) imaging is ~ 135 pc (Figure 
1), and even the most conservative estimates for the SFR re- 
sult in densities of w 3000 Mq yr^^ kpc~^. Fitting to the 
WISE 12 and 22/xm photometry, the integrated 8-1000 /xm lu- 
minosity of SDSS J1506-I-54 is estimated to be in the range 
\og{Li-[i/ Lq) = 12.6-13.1, assuming a representative range 
of spectral energy distributions appropriate for star-forming 
galaxies (Chary & Elbaz 2001, Dale & Helou 2002, Rieke et 
al. 2009). The dispersion in rest-frame 24fim luminosity us- 
ing this range of templates is just 0. 1 dex. A complementary 
estimate of Lm is therefore made using L24 using the em- 
pirical correlation between £24-^8- 1000 found by Rieke et 
al. (2008); we find log(LiR/L0) = 12.78 ± 0.06, consistent 
with the estimate above. The SFR of SDSS J1506H-54 is es- 
timated to be in the range 340- 1400 7\f0 yr^^ assuming the 
calibration of Kennicutt (1998), scaled to a Chabrier initial 
mass function. This is in good agreement with the SFR esti- 
mated from stellar population fits to the UV-optical-IR pho- 
tometry (DS12). 

The stellar mass is estimated from stellar population tem- 
plate fits to the 0. 1 -3 /^m photometry, with \oglQ{^I^, / Mq) — 
11.12 ± 0.07 (assuming a 0.1-100 Af© Chabrier initial mass 
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Figure 1. (left) Velocity integrated (AV = 360 kms^^) IRAM map of the CO(2-l) emission of SDSS J 1506+54. The white ellipse shows the 4.4" X 2.7" 
beam, and contours are at levels of 0.33 mjy (equivalent to Icr; negative values are shown as dashed contours), (right) Zoom-in of 5" X 5" box (left) showing 
the HST/WFC3 F814W (rest-frame V-band) image. The light is dominated by the compact core, re = 135 pc (DS12), but fainter, diffuse optical emission is 
also visible over scales of several kpc indicating a low surface brightness (LSB) envelope of UV emission (contours are shown as a guide). More extensive LSB 
emission beyond this field of view is also evident, at a level of /^f814w ~ 25 magarcsec^"^, including a potential tidal tail feature (P. H. Sell et al. in prep). 
This points to a merger origin for SDSS J1506-I-54 (§6). 



function, see Moustakas et al. 2013). Finally, the outflow ve- 
locity of the galactic wind is traced by the interstellar medium 
Mg 11 AA2796, 2803 absorption lines, which indicate a maxi- 
mum -i^abs ~ -1210kms-i (Tremonti et al. 2007, DS12). 

2.1. Comments on possible AGN contribution 

The source does not appear to be dominated by an 
AGN. Although SDSS J1506H-54 has the most luminous 
[Olll]A5007 and [Nev]A3426 emission of the DS12 sample 
(log(L[oiii]/ergs-i) = 42.1 and log(L[NoV]/ergs"i) = 
41.1), only 4 counts were detected with Chandra in the 2- 
lOkeV X-ray band, corresponding to a luminosity of Lx = 
2.9^2 9 ^ lO'''^ erg s^^. Assuming typical ratios between X- 
ray, [O lll]A5007 and [Ne v] A3426 luminosities, and account- 
ing for possible small ^30% dust attenuation of the narrow- 
line region based on the reddening derived from our fit to the 
SED (DS12), the Chandra limits imply an X-ray absorption 
factor of « 10-20 for a buried AGN (DS12, Heckman et al. 
2005, Gifli et al. 2010). However, as DS12 argue, even if the 
X-ray attenuation factor was as high as 100, the AGN would 
contribute less than half of the observed 22/im luminosity 
(Diamond-Stanic et al. 2009, Gandhi et al. 2009). Thus, star 
formation is likely to be the dominant power source of the 
bolometi-ic emission of SDSS J1506H-54. 

3. IRAM PLATEAU DE BURE 2 MM OBSERVATIONS 

Observations were conducted with IRAM PdBI over June- 
July 2012 in configuration D, using 5 antennas and the WideX 
correlator We targeted the CO(2-l) 230.54GHz rotational 
transition at I'obs = 143.37GHz (2mm band) in the direc- 
tion of SDSS J1506H-5402. The total (on source) integration 
time was 6.2 hours (after flagging; scans on-source were dis- 
carded for which the phases deviated by more than 45 degrees 
from the solution), during which time the system temperature 
ranged between T^ys = 100-400K ({Tsys) « 200-250K) 
for precipitable water vapour in the range pwv = 4-12 mm. 



Bandpass calibrators were the sources 3C84, 2200H-420 or 
3C273, gain calibration was performed with the sources 
1418H-546 and J1604H-572 and the source MWC349 was used 
for flux calibration. The flux calibration accuracy is ~5-10% 
at 2 mm. Data were calibrated, mapped and analyzed using 
GILDAS (Guifloteau & Lucas 2000). 

4. RESULTS 

Figure 1 shows the velocity integrated CO(2-l) map of 
SDSS J1506H-54 and Figure 2 presents the 2 mm spectrum 
around j^obs = 143.37GHz. The CO(2-l) line is detected 
with high confidence at >8(T in the integrated map. The 
integrated line flux is SAV = 0.97 ± 0.19Jykms~\ cor- 
responding to a luminosity of -^co(2-i) ~ ('^■^ -'- '-'■^) ^ 
10^ Kkms~^ pc^. Uncertainties are conservative; they are 
estimated via a bootstrap analysis, where we generate a se- 
ries of realizations of the spectrum, adding noise to each 
40km s^^ channel, selected from a Gaussian distribution with 
(7=1.1 mJy. The standard deviation of the ensemble of in- 
tegrated spectra are taken to be the 1 a error of the measured 
line flux. 

The line profile is adequately fit by a single Gaussian pro- 
file (x'^/iy = 1.1), with a velocity width FWHM = 283 ± 
31kms~^ (Fig. 2). Assuming the average inclination angle 
i — 30°, the true width could be as large as ^560 kms^^. 
Galaxy-integrated CO spectra often exhibit double peaked 
line profiles, indicative of a rotating disc or ring (e.g. Downes 
& Solomon 1998). There are some hints of this in the cur- 
rent spectrum, and so we also fit a double-Gaussian profile 
(fixing the individual line peaks and individual dispersions to 
be equal). The fit is formally x^/^ = 0-9> not significantly 
improved over the single component fit. The peak-to-peak 
velocity is AV = (147 ±12) km s^^, and the line dispersion, 
which includes the 1 -dimensional turbulent velocity disper- 
sion of the gas, is cr = (49 ± 6)kms^^. 

At first glance, it seems puzzling that the fine width is not 
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Figure 2. The 2 mm spectrum of SDSS J1506+54 obtained with IRAM PdBI, 
binned to a resolution of 40 km s~^. The channel noise (rm.s.) is l.lmJyper 
40 km s^ ^ channel. We consider two model fits to the data: a single Gaussian 
profile, and a double Gaussian (§4). The latter is an approximation to the line 
profile of a rotating disc or ring, and the observed spectrum does show hints 
of a double-peaked profile. We note that there is a 3.2(t spike in the spectrum 
close to the outflow velocity (as traced by blue-shifted Mg II absorption; we 
show the best-fit absorption profile of the Mg II doublet as a shaded region on 
an arbitrary flux scale). We acknowledge this is currently far too tentative to 
draw any conclusions, but we remark that - if real - this could be indicative 
of cold gas entrained in the outflow (§5.4). 

even larger, if the mass distribution is as compact as suggested 
by the rest-frame UV morphology. There are two main pos- 
sibilities: (i) the Te for the stellar mass is more widely dis- 
tributed than Te of light measured in the F814W band (near- 
infrared observations at the same resolution will be required 
to test this) such that the environment dominated by star for- 
mation is compact (r w re), but molecular gas dominated, or 
(ii) the CO-emitting gas is more widely distributed than the 
compact optical emission. The lack of size and orientation 
constraints for the CO(2-l) emission makes any estimate of 
A/dyn uncertain, however it would be unusual if the bulk of 
the gas reservoir (and presumably the concomitant dust) was 
beyond and configured in a way that did not completely 
obscure the UV-bright core. 

5. INTERPRETATION AND DISCUSSION 

5.1. Comparison with other galaxies 

The measured Lih/L'qq ss 1500 is one of the highest mea- 
sured for any star-forming galaxy (Fig. 3). Even if we ac- 
counted for an AGN contribution of 25% to our lowest es- 
timate of iiR, SDSS J1506H-54 is at the extreme tail of the 
Lik/L'qq distribution, with Lir/^qo ~ 550 in this case. 
There are a handful of other star-forming galaxies with similar 
Lik/L'qq properties; for example, the z = 0.575 hyper-LlRG 
IRAS F00235H-1024 (Combes et al. 2011) and the z = 0.633 
ULIRG IRAS F10398H-3247 (Combes et al. 2012) both have 
CO line luminosities measured in the 2-1 transition, and are 
also at the extreme end of the Lir/L^q distribution (Fig. 3). 
Several galaxies in the sample of Combes et al. (2012) have 
L'qq upper limits implying extreme Lir/X^q and of these, 
two are obviously hosting AGN, as evident from their optical 
spectra. 

Like SDSS J1506H-54, F00235H-1024 is not dominated by 
an AGN, although there is thought to be a non-negligible 
(w30%) AGN contribution to its total infrared emission 
(Verma et al. 2002, Farrah et al. 2002). F10398H-3247 could 
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Figures. Lir-Lqq (Jup < 3) plane for star-forming galaxies, showing 
the compilation of data for super star clusters, normal discs, local starburst 
galaxies, ULIRGs, submillimeter galaxies, hyper-LIRGs and circumnuclear 
starbursts from Andrews Sl Thompson (201 1). We plot L'qq Jup > 1 uncor- 
rected for sub-thermal excitation, but horizontal lines indicate how the points 
would move for corrections of 0.8 < r'21 < 1 and 0.5 < r^i < 1 re- 
spectively. The shaded region shows Lir «i ^Edd for optically thick dust 
emission (see §5.2). We also show predictions from Geach & Papadopou- 
los (2012) which models the global gas reservoir as a combination of dense 
star-forming gas and a diffuse quiescent cold gas phase. SDSSJI506+54 is 
close to the Eddington limit, with only a handful of other galaxies exhibiting 
similar ratios (Combes et al. 201 1 & 2012). 

also contain a deeply buried AGN, as it has a deep silicate 
absorption feature (ts, > 3, Dartois & Munoz-Caro 2007), 
often associated with nuclei blanketed by a dense screen of 
carbonaceous and silicate grains. The key difference between 
these galaxies and SDSS J1506-I-54 is in optical morphol- 
ogy: F00235H-1024 and F10398H-3247 are interacting sys- 
tems, where multiple UV/optical clumps can be identified in 
//5r optical and NIR imaging (Stanford et al. 2000, Farrah et 
al. 2002, Combes et al. 2012), thus it is unclear whether these 
galaxies support such extreme S*. 

5.2. Eddington-limited star formation 

How do we interpret such extreme systems? Can Xir / L'^q 
be driven to such high values through star formation alone? 
In individual galaxies, isolated star-forming cores might be 
Eddington limited, which is to say is capped by radiation 
pressure from the recently formed O and B associations (e.g. 
Scoville 2004, Shirley et al. 2003). However, when taken as 
a whole, galaxies form stars at sub-Eddington rates due to in- 
termittency - the fact that only a small fraction of the cold gas 
reservoir is actually forming stars (Downes & Solomon 1998, 
Murray, Quataert & Thompson 2005, Andrews & Thompson 
201 1). Under certain conditions however, one could envision 
a scenario where the majority of the gas reservoir could be 
driven to form stai-s at the Eddington Hmit. SDSS J1506H-54 
is an excellent candidate for this phenomenon. 

In the absence of nuclear heating, and assuming the model 
of momentum-driven feedback (e.g. Murray, Quataert & 
Thompson 2005) the upper limit of Xir/ L'qq for star-forming 
galaxies is set by the Eddington luminosity. For optically 
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thick dust emission (rioo/^m > 1): 

^Edd = XcoL'qq (1) 

K 

(Andrews & Thompson 2011), where here we set Xqo = 3 
(see §5.3) and k is the Rosseland-mean dust opacity. In Fig- 
ure 3 we show the Eddington luminosity predicted for k = 5- 
lOcm^g-i iUs/f^^), where f^^ = i/iso is the Milky 
Way dust-to-gas mass ratio (see Andrews & Thompson 201 1). 
We assume a dust-to-gas ratio of /dg ~ ^/50, appropriate for 
dusty star-forming galaxies (e.g. Kovacs et al. 2006). Increas- 
ing either /dg (as might be appropriate in highly obscured 
galaxies), or k, or decreasing Xqo, lowers Z^Edd for a fixed 
^'co- 

Our chosen values aim to be conservative in this respect, 
such that no star-forming galaxies can be described as super- 
Eddington, but we caution the reader that in this parame- 
ter space the Eddington limit is not expected to be a hard 
edge. A simple framework that explains the intrinsic range 
in Lir/Lqq is the model of Geach & Papadopoulos (2012) 
that predicts the minimal dense gas mass for a given Lir , and 
a range of gas mass ratios ^sf = Mdcnsc/Mtotai- Typically, 
^gp « 0.02 for quiescent discs and ^sf ~ 0.2 for ULIRG- 
like systems (Fig. 3). The Eddington limit in equation (1) is 
approached as fsF 1, indicating that the majority of the 
gas in SDSS J1506H-54 is in a dense, active phase and the in- 
termittency is close to unity. 

5.3. Estimating the gas mass 

Recent studies of the local (U)LIRG population where well- 
sampled CO and ^■^CO ladders are available (Papadopoulos et 
al. 2012a) show that, when high-/ CO or heavy rotor (HCN, 
CS) line transitions are included in radiative transfer mod- 
els, the A/h2 / L'qq = Xco factor increases by a factor ^3- 
10 compared to the widely applied Xqo = 0.8 conversiorQ 
usually adopted for ULlRG-like systems (Papadopoulos et al. 
2012b). The explanation is that in the supersonic turbulent 
discs of ULIRGs, a significant fraction of the gas mass is ex- 
pected to be in a dense (n > 10'* cm^'^) phase (Padoan & 
Nordlund 2002, Papadopoulos et al. 2012, Papadopoulos & 
Geach 2012). This is not well traced by single low-7CO trac- 
ers, despite the fact that these were originally used to cali- 
brate the value of Mh^/L'^^q in ULIRGs (Downes, Solomon 
& Radford 1993, Solomon et al. 1997, Downes & Solomon 
1998). 

Here we assume that the majority of the gas reservoir is 
indeed in a dense, turbulent phase. Setting Xqo — 3 and 
r2i — 0.8 (Papadopoulos et al. 2012a, Geach & Papadopou- 
los 2012), we estimate Mu^ = (1.9 ± 0.3) x IO'^^Mq. This 
is in agreement with the minimum dense gas mass expected 
if LiR is powered by Eddington-limited star formation, since 
one expects a maximum e^, = im/Afdcnsc ~ 500 {Lq/Mq) 
(Scoville 2004, Thompson et al. 2005, Thompson 2009). The 
e,t expected in feedback models is indeed close to the value 
actually measured both for resolved star-forming cores in spi- 
ral discs and even entire starbursts (Scoville 2004, Shirley 
et al. 2003). Similar values of e,t are obtained when one 
considers the mass-to-light ratio of a deeply dust-enshrouded 
(tir > 1) zero-age main sequence (ZAMS) population, with 
-^ZAMS ~ ^IR (Downes & Solomon 1998). It is important 

^ we omit the units of A/0 (K km s^^pc^)^^ 



to note that, regardless of the exact gas mass calibration, as- 
suming that the majority of the infrared emission originates 
from enshrouded star-forming regions, the empirical scalings 
in Figure 3 clearly indicates that SDSS J1506H-54 is close to 
the Eddington limit. 

5.4. Hints of gas entrained in a wind 

The compact starburst in SDSS J1506H-54 should be effec- 
tive in the launching of a powerful momentum-driven galactic 
wind, powered by radiation pressure on dust grains intermin- 
gled with the ISM, and the ram-pressure of supemovae deto- 
nations (Murray, Quataert & Thompson 2005). The strongly 
blue-shifted Mg II lines are excellent evidence that this is the 
case (Tremonti et al. 2007, DS12). 

The 2 mm spectrum of SDSS J1506H-54 exhibits a 
marginally significant (3.2(t) narrow (cr < 40kms~^) spike 
close to the outflow velocity traced by Mg II (no RFl or tel- 
luric contaminant is expected at this frequency). We can draw 
no conclusions on such a low-significance feature, but we re- 
mark that if cold gas has become swept-up in the outflow 
traced by Mg II absorption, and have survived shock-heating, 
then we might expect to detect entrained cold gas in this way. 

6. SDSS J1506+54 AS AN EDDINGTON-LIMITED STARBURST 
INDUCED IN A LATE-STAGE MERGER 

6.1. Formation 

We speculate that SDSS J1506H-54 is the remnant of a major 
merger In this picture angular momentum losses caused by 
the dissipation of energy in shocks and gravitational torques 
cause the collapse of disc-gas to the central regions in each 
parent, and subsequently into a single merged core. This is a 
classic picture of galaxy growth that has been supported for 
many years by both theory, numerical simulations and obser- 
vation (Toomre & Toomre 1972, Sanders et al. 1988, Hern- 
quist 1989, Barnes & Hernquist 1991, Sanders & Mirabel 
1996, Schweizer et al. 1996, Mihos & Hernquist 1996, Wuyts 
et al. 2010, Hopkins et al. 2012). 

Tidal interactions in earlier stages of the merger might well 
be responsible for the faint extended (over several kpc) op- 
tical light (Figure 1, P. H. Sell et al. in prep). Episodes of 
massive star formation must have occurred in the recent past 
(0.5-1 Gyr), with the A stars formed in (or moved to) rela- 
tively unobscured regions, in order to produce strong H(5 ab- 
sorption observed in the optical spectrum. Indeed, the pres- 
ence of post-starburst spectral features is entirely consistent 
with the typical timescales for equal-mass mergers (Lotz et 
al. 2008). 

6.2. Fate 

In the absence of additional cooling of gas, SDSS J1506H-54 
must be close to the cessation of star formation. At the cur- 
rent rate the gas will be consumed within a few tens of Myr, 
increasing the stellar mass by just ^15%. This will place the 
galaxy in the exponential tail of the stellar mass function at 
z K, 0.5 (Moustakas et al. 2013), but the current starburst 
contributes only a modest fraction of the total mass. Never- 
theless, this small addition can have important consequences 
for the properties of the descendant (Robaina et al. 2009, Hop- 
kins et al. 2010, 2012, Hopkins & Hernquist 2010). 

Signatures of compact Eddington-limited starburst episodes 
could be found in the 'fossil record' of local spheroids. Inte- 
gral field observations of the cores of local elliptical galax- 
ies reveal examples of compact (^^lOOpc), 'kinematically de- 
coupled', relatively young ( ^ 5 Gyr) stellar components in 
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a subset of the population (e.g. McDermid et al. 2006). The 
origin of these structures could well be found in powerful, 
compact starbursts such as the one we present. Star formation 
driven at extremely high could also have a critical impact 
on the form of the stellar initial mass function, due to the cor- 
responding high cosmic ray densities in such environments 
(Papadopoulos et al. 201 1). 

7. SUMMARY 

We have presented observational evidence of what we term 
a 'redline' starburst: a galaxy forming stars close to the 
Eddington limit assuming current models of stellar feed- 
back. The measured Lir/L^q ^ 1500 is one of the high- 
est measured for any galaxy. Our conclusion is that the 
SDSS J1504H-54 system is a final stage merger, undergoing 
a high intensity circumnuclear starburst. If we are witnessing 
the intense - but fleeting - closing stages of the assembly of 
a massive stellar bulge, SDSS J1506H-54 represents a unique 
opportunity to study the physics associated with this impor- 
tant evolutionary phase, and the mechanics of star formation 
at its most extreme. 
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